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Abstract 

Within the mode-coupling theory for structural relaxation in simple systems 
the asymptotic laws and their leading-asymptotic correction formulas are de- 
rived for the motion of a tagged particle near a glass-transition singularity. 
These analytic results are compared with numerical ones of the equations of 
motion evaluated for a tagged hard sphere moving in a hard-sphere system. 
It is found that the long-time part of the two-step relaxation process for the 
mean-squared displacement can be characterized by the a-relaxation-scaling 
law and von Schweidler's power-law decay while the critical-decay regime is 
dominated by the corrections to the leading power-law behavior. For param- 
eters of interest for the interpretations of experimental data, the corrections 
to the leading asymptotic laws for the non-Gaussian parameter are found to 
be so large that the leading asymptotic results are altered qualitatively by the 
corrections. Results for the non-Gaussian parameter are shown to follow qual- 
itatively the findings reported in the molecular-dynamics-simulations work 
by Kob and Andersen [Phys. Rev. E 51, 4626 (1995)]. 

PACS numbers: 64.70 Pf, 61.20 Lc 
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I. INTRODUCTION 



The mode-coupling theory (MCT) for the structural relaxation of glassy liquids |],|] has 
provided explanations for long-known phenomena as for example the stretching of the a pro- 
cess, the time-temperature superposition principle, and the von Schweidler law. It has also 
predicted new phenomena such as a square-root singularity in the temperature dependence 
of the Debye- Waller factor and the critical decay. For details the reader is referred to Ref. 

and the papers quoted there. The universal features of the MCT arise as leading-order 
results near a bifurcation singularity, which was identified as an ideal glass transition — in 
reality smeared out due to additional relaxation processes. These theoretical findings also 
provided some motivation to a number of experiments focusing on the evolution of structural 
relaxation in the vicinity of the predicted glass-transition singularity (see Ref. P-IPTf and 
the papers quoted there). In turn these experiments have sparked off new theoretical devel- 
opments such as the derivation and discussion of next-to-leading-order asymptotic results 
in the vicinity of the bifurcation singularity |18[ and the extension of the theory to liquids 
consisting of or containing non-spherical molecules |T9| , pO| . Using the coherent intermediate 



scattering function of a hard-sphere system (HSS) as an example, the preceding work |18 



illustrates how asymptotic corrections explain the deviation from the leading-order results 
and shows how general properties of these corrections may be included into a quantitative 
analysis of experimental data. Due to the presence of additional amplitudes and a new 
second-order scaling function, many features of these corrections are different for different 
observables, used for the description of the glassy dynamics. Therefore we intend to ex- 
tend the previous work to a discussion of the conceptually simplest functions characterizing 
the dynamics, namely the density correlator of a tagged particle and its special limits, the 
mean-squared displacement and the non-Gaussian parameter. 

The tagged-particle density correlators can be measured by a variety of techniques. Neu- 
tron scattering from molecular liquids with an incoherent scattering cross-section such as or- 
thoterphenyl fP7| , |2"T|1 probes its spectrum. For a glassy colloidal suspension the correlator has 
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recently been measured by dynamic light scattering . From measurements of the incoher- 
ent intermediate scattering function for small wave vectors one can extract the mean-squared 
displacement P2-2^| and, in principle, also information about the non-Gaussian parameter. 



In computer simulations the tagged-particle correlator is a preferred quantity — compared to 
the coherent density fluctuation function — due to statistical advantages. The mean-squared 
displacement and the non-Gaussian parameter can also be directly extracted from computer 



simulation data [ njPql ■ The reported simulation results for the non-Gaussian parameter 
challenged the universal features predicted by MCT. Within the accessible simulation-time 
window the non-Gaussian parameter showed no signs of the predictions of a two-step re- 
laxation process nor of the time-temperature superposition principle, even though other 
dynamical quantities such as the intermediate scattering function, self-intermediate scatter- 
ing function and the mean-squared displacement fitted into the picture for glassy dynamics 
drawn by the MCT. 

To proceed, equations of motion will be derived for the mean-squared displacement and 
the non-Gaussian parameter. The relaxation kernels in those equations are given by mode- 
coupling functionals, which require the correlators of the density and the tagged-particle 
density as input. For the latter the previously-derived MCT equations are used. Then 
an asymptotic expansion of the solutions will be carried out. All equations will be solved 
and studied in quantitative detail for the model system of a tagged hard sphere immersed 
in a HSS. The HSS has already been the subject of a number of theoretical investigations 
as can be inferred from Refs. |T8| , |26[| and the papers quoted there. Its importance stems 
in part from the fact that it is the simplest system for which a glass transition has been 
detected experimentally |]8|j9|j27[| , thus providing the archetype for quantitative tests of the 
theory. Comprehensive comparisons between MCT results for the evolution of structural 
relaxation and the corresponding experimental findings obtained by dynamic light scattering 
for hard sphere colloids have been published in Refs. |f|j9[ p8| , p5| .In the data analysis enter 
two numbers as fit parameters: a time scale t , given by the viscosity of the solvent, and the 
critical packing fraction ip c . The predicted value for ip c || differs from the the experimental 



value by about 12%. Here as in any other test of a singularity theory, data have been studied 
as a function of the distance e — (tp — <p c ) / (p c from the position tp c of the singularity. The 
cited tests have been summarized as follows |30|J3TJ| : the leading asymptotic MCT results 
account for the data within the experimental uncertainties for the liquid states (e < 0) and 
also for the glass states (e > 0). This does not only hold for universal formulae like the 
scaling laws. It holds also for the numerical values predicted for the exponents, the master 
functions entering the scaling laws, and for the various wave-vector-dependent amplitudes. 
Therefore it seems justified to illustrate our theory for the tagged particle motion in a 
HSS. The cited and tested results for the density dynamics will enter as input taken from 
the preceding work |fL8|j . But the HSS is also relevant because its dynamics shares many 
features with other simple systems as will again be demonstrated in this article. 

The paper is arranged as follows. After briefly reviewing the basics of the MCT glass 
transition (Sect. H) we start by discussing the properties of the tagged-particle density 
correlator in Sect. [ITT|. The next Section covers the mean-squared displacement and its 
asymptotic behavior. The quality of the Gaussian approximation for the self-intermediate 
scattering function of the HSS is assessed (Sect. [V]) before turning to the non-Gaussian 
parameter in Sect. [VT|. In the conclusion we will summarize our results and suggest how 
some recent experiments and simulations on single-particle quantities can be interpreted. 



II. THE IDEAL MCT GLASS TRANSITION 

The idealized MCT deals with the structural dynamics of simple liquids, consisting of N 
particles at positions fj, j = 1, . . . , N. A self-consistent treatment of the cage effect, which 
is thought to be the reason for the glass transition, leads to a closed set of integro-differential 
equations for the density correlator $ g (i) = (p^(O)p^(t)) / S q , i.e., the autocorrelation func- 
tion of the density fluctuation Pq{t) = 1 fy/N ^ . exp [iqfj (t)] for wave vector q, normalized 
by the static structure factor S q = (p|-(0)pq(0)); () signifies an average with respect to the 
canonical ensemble and q = \q\ abbreviates a wave vector modulus. Specializing to a colloid 
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model, the MCT equations of motion read 

t 

r q % (t) + % (t) + Jm g {t- t') $ g (*') dt' = 0. (1) 
o 

They are specified by the times r q = S q / (D q 2 ) with D denoting the single particle diffusion 
coefficient. The memory kernel in this equation is given as a mode-coupling functional 
m q (t) = T q (<3> (t)), where T q is determined by the structure factor S q [Q], which depends 
smoothly on external control parameters such as the density n. 

At a critical density n c the above equations exhibit a bifurcation in the long-time limits 
f q = lim^oo $ g (t) of the density correlators. The bifurcation singularity can be identified 
with a liquid-glass transition: At n c , the Debye-Waller factor f q , also called form factor, 
jumps from zero to the critical from factor fE > 0. The instability of the glass for n — > n c + is 



reflected by a square-root singularity of the Debye-Waller factor: f q = f q +h q ^Ja/ (1 — A) + 
0{a). Here, a = Ce is called the separation parameter, 0.5 < A < 1 is called exponent 
parameter, and h q = (l — f q ) 2 e q > denotes the critical amplitude. The quantities e q , 
A, and the constant C > 0, relating a and the reduced density e = (n — n c ) /n c , can be 
calculated by straight-forward formulae from the functional T 0. 

Via a simple transcendental equation the exponent parameter A determines two anoma- 
lous exponents: the critical exponent a, < a < 0.5, and the von Schweidler exponent b, 
< b < 1. Furthermore A fixes a constant B > 0. There appears a single time to, specifying 
the scale for the transient dynamics. The bifurcation dynamics is then ruled by two critical 
time scales, denoted by t a and t' a : 



t ^ io/ i^, c = toB -v Vk r,, = i-, 7 = i- + i. (2) 

The mathematical relevance of these concepts is evident from the following limit results || . 
First, lim t -^ lim^o j [®q (ito) — f q ] / t a = hq- Thus, near the transition the plateau 
value f q is approached from above according to a power law, called critical decay 

% (t) - H ~ h q {t /t) a , t < t < t a . (3) 



Second, 

lim^olim^o- [* q (K) -ff\/t b = ~K = lim^ 00 lini ff ^o-{[^(tt ff )-/, c ]/vW}/* & - 
Near the transition the plateau value /° is left in the liquid state according to a power 

law, called von Schweidler decay 

% (t) - f c q - -K (t/t'J , t a « t « C (4) 

There are two dynamical scaling laws describing the bifurcation dynamics near the tran- 
sition in a leading-order asymptotic limit. The first one describes the dynamics on a scale 
t a : $ q (t) - j c q ~ h q G (t). Here G (t) = (t /t) a for a = and 

G(t) = yR<7±(t/t CT ),<7^0. (5) 

The a-independent master functions g± are determined by A. They describe the crossover 



from the critical decay g±(i <^1) = l/i a to arrest in the glass, g + (£ 3> l) = \/y/l — A, 
or to von Schweidler's decay in the liquid, g_ (t ^> l) = — -Bt 6 + (Bt b y The result 
is obtained by solving the MCT equations using |$ 9 (t) ~ fq\ as small parameter. This is 
equivalent to writing t = it a and then expanding in powers of the small parameter y|£t|. 
Pushing the expansion to next-to-leading order extends the first scaling-law result to 



% (t) = f q + h q G (t) + h q \H (t) + K q G (ty + aK q ^ . (6) 

Here two further amplitudes K q and K q and a new function H (t) appear. The latter obeys 
a scaling law: H (t) = k (a) (to/t) 2a for a = and H (t) = \a\ h± (t/t a ) for o ^ 0. All these 
quantities can be evaluated from the mode-coupling functional T q as is comprehensively 



explained in Ref. |18|], where we changed notation Kto K q . The domain of applicability of the 
expansion (||) is called /3-relaxation window, and correspondingly G and t a are respectively 
referred to as f3 correlator and /5-relaxation time scale. 

The second scaling law describes the liquid dynamics on scale t' a . It deals with the decay 
of the correlator from the plateau to zero and it is based on the relation lim CT ^o- &q y^a) = 
$q (t) . Here the a-independent functions $ g (t) are the solution of the implicit functional 
equations 
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J [m, (If) - (*')] ^ = y m, (* - $ 9 (*') (7) 



o o 



The kernel fh q (tj = T q (t)j is the mode-coupling functional at the transition point 
and the equation is to be solved with von Schweidler's law as initial condition. Near the 
transition within the window t a <C t this implies in leading order 

% (t) = % (t/Q . (8) 

This window is referred to as the a regime, t' a is called a-relaxation time scale, and Eq. (||) 
is the superposition principle for the a process. 

Every theory dealing with the dynamics of a variable A as a probe of structural relax- 
ation in glassy systems relies on an understanding of the dynamics of density rearrangements. 
Hence a microscopic theory will require results for $ g (t) as input. For the following dis- 
cussion it is therefore necessary to appreciate the concepts and results formulated above. 



They are explained and illustrated in Ref. |TB| for the HSS. The quantitative examples of the 



following discussion shall also be done for that model. We follow the previous conventions 
by using the particle diameter d as unit of length, d — 1, choosing the time unit so that the 
short time diffusivity Dq = 1/160, and using the packing fraction (p = 7md 3 /6 as the control 
parameter. In this case one gets: tp c = 0.516, C = 1.54, A = 0.735, a = 0.312, b = 0.583, 
B = 0.836, B x = 0.431 ,t = 0.425. 

III. THE TAGGED PARTICLE DENSITY CORRELATOR 

A. Equation of motion 

The dynamics of a tagged particle is explored by its density correlator (t) = 
(pp (0) P%(t))i where p%(t) = exp [iqr (t)]; r denotes the position of the particle. For a 

system driven by Brownian dynamics the MCT equation of motion reads 

t 

« (t) + $ s q (t) +Jm s q (t- & q (f) dt' = (9) 
o 



with r* = 1/ (D^q 2 ). Here _Dg * s the short-time diffusion coefficient of the immersed parti- 
cle. The memory kernel m q (t) = T q (<& if) , $ s (i)) is expressed through the mode-coupling 
functional J 7 ^, where not only (t) but also $ g (£) enter: 

^ (/, f s ^=jtfS nSkcf (?) A (io) 

Here = (p|?(0)pg-(0)) / (nS^) denotes the single-particle direct correlation function 0. 
To solve Eqs. @ and (|T0|) numerically we proceed as in Ref. | 18fl : The mode-coupling 



functional is rewritten in bipolar coordinates by change of variable and the two remaining 
integrals are approximated by Riemann sums 

r q (f, n = nA 3 / (le^vr 2 ) j2 E ' 3 >> (^/^ 5 ) ( p + ^ 2 - 2 ^/*/;- ( n ) 

The prime indicates that the sum runs only over those values p for which k, p, and q obey the 
triangle inequality. Thereby, we again get a precisely defined tagged-hard-sphere-particle 
model described by 100 coupled integro-differential equations for (t) on a grid of 100 
equally spaced wave numbers q = Aq with q = ~, . . . , 99^ and step size A = 0.4. 

The substitution of Eq. (|TT| ) for the mode-coupling functional (|10D might seem rather 
crude. However we have checked that the results seldom differ by more than a few percent 
from those obtained with 300 or 900 grid points for the Riemann sum. Only at small wave 
numbers qualitative changes occur due to bad angle resolution; but these do not influence 
other wave vectors because the major contributions to the integral in Eq. (|1(]) come from 
the wave vectors near the structure-factor-peak position. 

In our application to the HSS the static structure factor S q and the direct correlation 
function c s q are calculated in the Percus-Yevick approximation [^]. Assuming short-time 
diffusion according to Stokes' law, which entails Dq/Dq = d/d s , only two control parameters 
remain: the packing fraction ip of the host particles and the diameter d s of the tagged 
particle. All results for the coherent density correlators $ g (t) of the HSS, which will be 
needed in the following, will be taken from the solution of Eq. ([!]) discussed in Ref. []18 . 
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B. The critical Lamb Mossbauer factor 



The trapping of a tagged-particle by its surrounding host particles manifests itself by a 
nonvanishing Lamb-Mossbauer factor /* = lim^oo (t) . From Eq. @ one finds that 
/* is a solution of 



^ (/,/*)■ (12) 



It is distinguished from other solutions of Eq. flT2|), say /* by the maximum property, /f > 
fg ||. It can be found by the iteration: /* = ]xni n -> 00 f^ n \ where /q^ n+1 V ^1 — fq^ n+1 ^ = 
(/; f s ^)i fq^ — 1 HI]- Obviously / = implies / s = 0, but not vice versa. The 
spatial Fourier transform f s (r) of the Lamb-Mossbauer factor is the probability of finding 
the particle for t = oo at distance r from where it started for t = 0. Therefore, a particle 
cannot be trapped as long as the host particles are in a liquid state. But even when the host 
particles are in a glass state the tagged particle may still be able to diffuse in the arrested 
structure. The various scenarios arising from Eqs. (|12D and (0) have been discussed before 



(compare Ref. |34| and the papers cited there). 

For a tagged hard sphere of diameter d s = 1.0 the Lamb-Mossbauer factor jumps from 
to the critical value f^ c > 0, which is shown in Fig. f[[ when the packing fraction (p passes 
through the glass-transition singularity ip c . For d s = 0.6 the particle still gets trapped at ip c 
but the /| c -versus-g graph is narrower than before (see Fig.^j) . This means that the centre 
of the smaller particle can explore a larger volume. Reducing the particle size d s further, 
the system passes through a percolation threshold at some critical diameter d sc . In the 
following we will only be interested in tagged particles with diameters d s > d sc , for which 
the Lamb-Mossbauer factor shows the generic fold bifurcation. 



C. Asymptotic laws 
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1. P relaxation 



Close to the bifurcation the correlators <3>* (t) show a characteristic two-step relaxation 
process, as can be seen in Fig. |3|. In Ref. || it was discussed and in Ref. [35) it was 



demonstrated for relevant examples that the long-time behavior of correlators, which are 
the solutions of MCT equations such as (ffl) and (|^), is determined by 



1 - s$°(s) 

Here $*(s) denotes the Laplace transform $ (s) = C [$ (t)] (s) := / °° e~ s *$ (£) of (£). 
These equations depend only on equilibrium quantities via the static structure factor. They 
are therefore independent of the microscopic dynamics and any additional regular contribu- 
tions to the kernels m q (t) and m q (t) . Consequently, it makes no difference for the asymp- 
totic results discussed throughout this paper that we chose Brownian instead of Newtonian 
dynamics and omitted regular contributions to the kernels. 

Inserting the asymptotic expansion @ into Eq. (pT3[), using Taylor's theorem on the 
mode-coupling functional T ql and collecting like-order terms one obtains in leading and 
next-to-leading order 

K (*) = /r + K G (*) + K \ H (t) + K « G (t) 2 + ■ (14) 



Here the critical amplitude h s = (l — f q c ) e| for the tagged particle correlator and the 
correction amplitudes K s q and K q are to be calculated from 

Efo*- c $ e * = E c £ e *> ( 15a ) 

k k 

E (v - C S) ^ = - (i - /r ) < x + E c £ e * K * 
k k 

+ E G qpk e p e k + E G> q,pk e P e k + E G qp,k e p e k, (15b) 
p,fc p,fc 

E - q c ) cj^e = - (i - /r ) < + ^q/c + e cr* e *^*- ( 15c ) 

Taylor coefficients of the mode-coupling functional are introduced by 
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c s = c sc t5r s f a- o (f 2 \ 

qp\...p n ,k\...km qpi—Pn,ki...k m ~ qpi —PnM ...km ~ \ ) 

1 l d n d m F*(f c ,f sc ) 
n\ m\ df Pl . . . df Pn df kl . . . df km 

(1 " f£f •••(!- fZY t 1 - fkf ■ ■ ■ (1 - fkS ■ (16) 

For the derivation of these general formulae we have made no use of the linearity of the 
mode-coupling functional in / and f s . Figures [I] and ^] exhibit the HSS results for the 
amplitudes for the two diameters d s = 1.0 and d s = 0.6, respectively. Following Ref. ]T8] we 
show 



K s q = K s q VT^\ + Kl/VT^ -n + h + (oo) VT^ (17) 

instead of K^, where k = 0.961 and h + (oo) = —2.21 for the HSS. 

Equation (|14D for the tagged-particle correlator is formally identical to the above expan- 
sion (|D for the coherent density correlator. This reflects the universality of the asymptotic 
features of the MCT. Within the /5-relaxation window the dynamics is fixed by the f3 cor- 
relator G (t) in leading order and by G and H in next-to-leading order in an asymptotic 
expansion. Correlators referring to different probing variables merely differ in the various 
amplitudes such as (fq,fq C ), (h q ,hq), etc. These amplitudes do not depend on time nor on 
control parameters. However, there are rather important differences for the procedure to 
determine the amplitudes reflecting the fact that the bifurcation in the tagged-particle dy- 
namics is brought about by the bifurcation in the host-particle dynamics. The amplitudes 
e q , K q and K q for the cage-forming particles are the solutions of singular matrix equations 
and the solubility conditions in second and third order of the asymptotic expansion deter- 
mine the scaling functions G and H up to a single time scale to, which has to be matched 
at the critical point fL8f . By contrast, the positive matrix C q % has a spectral radius smaller 
than unity |33|]. Hence, the matrix S q k — C q ^ is regular and the amplitudes are uniquely 
determined by solving the linear equations (p!5|). The scaling functions G and H are solely 
determined by the properties of the host particles. The interaction potential of the tagged 
particle enters the Lamb-Mossbauer factor and the amplitudes via the direct correlation 
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function c s q , which is hidden in the Taylor coefficients (p!6[). No further time scale needs to 
be fixed. Even a change of the microscopic time scales r q in Eq. (|9[) would not affect the 
time scale of the asymptotic behavior. One can always go so close to the critical point that 
Eq. fll^D holds within a prescribed error margin. 

The first two terms on the right-hand side of Eq. are known as the factorization 

theorem: If the density correlators are rescaled accordingly the curves collapse onto the f3 
correlator G in an intermediate time window, the so-called /3-relaxation region, as is demon- 
strated in Fig. § for four wave numbers. Obviously, the /3-relaxation region depends on the 
wave number (upper panel) or — put more generally — on the quantity under consideration 
(lower panel). This is accounted for by the two g-dependent amplitudes K q and K q occur- 
ring in the leading correction to the factorization property. The (3 window expands as the 
system moves closer to the critical point as is demonstrated by a comparison of § s (t) for 
q<i = 10.6 for n = 9 and n = 10 in Fig. f|. More can be said if we consider that the dom- 
inant corrections are the short-time corrections to the critical law (|3|), and the long-time 
corrections to the von Schweidler law (f|) or the corrections to the Lamb-Mossbauer factor 
f q on the liquid (e < 0) or glass (e > 0) side, respectively. The short-time expansion of the 
right-hand side of Eq. (|TJ) is given by f s q c + h s q {t/t o y a {1+ [k (o) + K s q ] {t/t o y a } so that 
the factor k (a) + K q determines the sign and the strength of the leading correction to the 



critical decay. For the HSS one gets k (a) = —0.002 |18| and K q < K s v for q < p (compare 
Fig. |l]). This explains in particular why at early times the curves with smaller wave number 
fall below the curves with larger wave number. Since k (a) + K q changes sign between q2 
and <?3, the curves 0, 1, and 2 fall below and the curve 3 is above the dashed asymptote for 
t < 10. On the liquid side the analogous formula $^ (t) = f q c - h s q i b {l — [« (-b) + K*\ ?} 
provides the explanation for the deviations from the von Schweidler asymptote, Eq. (^), at 
later times: the g-dependence of the deviations is again determined by the amplitude K q . 
This entails that if one curve lies above the other when deviating from the factorization 
theorem at short times it will also do so at long times. This rule, which has already been 



discussed for $ ? (t) in Ref. Hl8| , holds for all quantities that asymptotically obey Eq. (|14|). 
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For the HSS one gets k (—6) = 0.569 [18| and therefore Fig. [I] implies that k (—6) + K q 
changes sign between qi and qi- This explains why the curves and 1 are below, and 2 and 
3 above the dashed asymptote in Fig. |] for log 10 t ~ 6. On the glass side the long-time 
behavior of § a q (t) is described by /* = f s q c + h s q y/a/ (1 - A) [l + yjd (R s q + «)] . 

Armed with Fig. ^| and what was explained for d s = 1, it is easy to infer what Fig. |] 
would look like for a hard-sphere of diameter d s = 0.6. Since it was already explained in 



Ref. JL8J how the deviations from the leading-order translate into the frequency domain, 
spectra shall not be discussed in this paper. 

In order to discuss the quality of the asymptotic expansions in a controllable manner, 
some convention for the error margin has to be made. There are various possibilities |18 



In this paper the range of validity shall be defined as the time window where for a given 
value $ of the correlator, the time t$ for the solution, $ = $(t$), deviates from the time 
t| s for the asymptotic expansion, say $ = f q + h q G(t$), by less than e per cent. This 
definition can always be applied to monotonic functions of time, in particular to functions 
in the /9-relaxation window. If one chooses e so small that the deviations from the leading 
order are quantitatively explained by the leading corrections, the range of validity for the 
leading-order /3-relaxation asymptote can be shown to expand with the (3 time scale t a for 
o — > 0+ and with the a time scale t' a for o — > 0—. In Fig. f| and some of the following 
figures the endpoints of the time intervals with e = 20% are marked by various symbols. 



2. a relaxation 

Within the a-relaxation window the tagged-particle correlator & q (t) obeys the super- 
position principle as formulated in Eq. (||) for the density correlator: (t) = <& s {t/t' a ). 
The control-parameter-independent master function (Pj is obtained from Eq. (|7|), where 
merely superscripts s have to be added to correlators and kernels ||. 

In Fig. |5| various rescaled correlators $^ (t) are compared to the a master function. The 
superposition principle best describes the a relaxation for later times. As the glass transition 
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is approached from the liquid side, the range of validity of the superposition principle, which 
is marked by diamonds and defined in analogy to the /3-relaxation regime above, extends to 
earlier rescaled times t. Close to the critical point, the deviations from the a master curve 
can be understood in terms of the leading-order correction 5<& s q (t) to the a-scaling law, 
which is linear in a: $j(t) = (?) +5& q (i)+0{a 2 ) 0. Although the equation governing 
the time evolution of 5<& s q (t) is rather involved, its short-time behavior can be expressed in 
terms of the amplitudes introduced above in connection with the (3 relaxation. One gets up 
to errors of order at b : 

5& q (?) = -ah\ {B 1 t- b + [k (-6) - 2B 1 K° - K s q ] } (18) 



with K,(—b) = 2.97 for the HSS [18|. The t~ b term, which diverges for t — > 0, dominates 

the deviations from a scaling for small \a\. Consequently, asymptotically close to ip c , the 

deviations occur at short times, where the a master function is described by von Schwei- 

dler's law (^). Therefore, the t-window for the a scaling expands to smaller rescaled times 
— 

proportional to |cr| 2& . For small |er| the a-relaxation window becomes independent of the 
correction amplitudes K q and K q and thus independent of the observable under consider- 
ation, because the observable-dependent amplitude h q occurs as a prefactor to both von 
Schweidler's law and the relevant correction t~ b . The dotted line in Fig. [| illustrates how 
the deviations 5<t q (tj from the a master function can be described by the leading correction 
-ah q B x t- b . 

IV. THE MEAN-SQUARED DISPLACEMENT 

A. Equation of motion 

The equation of motion for the mean-squared displacement 5r 2 (t) = (\r(t) — r*(0)| 2 ) 
36| can be obtained from Eq. (|9]) by exploiting its relation to the small-wave-number 



behavior of the tagged-particle density correlator (t) = 1 — q 2 5r 2 (t) /6 + O (q A )'- 
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z 

5r 2 (t) + D oJ m{0) (t - t') Sr 2 (t') dt' = 6D s t. (19) 



o 



Here we have introduced the kernel (t) = lim^o q 2r m s g (t)- Carrying out the limit q — > 
in Eq. (|T0|) one finds the representation of the new kernel as a new mode-coupling functional, 

m (o) (t) = T MSD ($(t), $ s (t)), where 



DC 



•^msd (/, /') = ^2 / nS k c s2 k*f k f k dk. (20) 



o 



Again this integral shall be rewritten as a Riemann sum over the previously introduced 
wave- vector grid of 100 terms 

99§ 

Tmsd (/, /') = nA 5 / (6tt 2 ) ^ Sktffrfkfi. (21) 

Two side remarks concerning the preceding formulae might be of interest. First, instead 
of solving Eq. (|1^) with <& q (t) and as an input for the determination of the kernel 

rrS°\ one could obtain the desired result directly as the small-wave-vector limit: Sr 2 (t) = 
6 linig^o [l — $o (t)] /q 2 - We have studied this procedure and checked the result to be the 
same as obtained from Eqs. ( |T9"D and (^). However, to do this, we had to use more than 100 
grid points for the Riemann sum and had to extrapolate carefully from small q to zero. This 
is necessary, since the discretization in Eq. ([TT| ) is too crude to produce reliable results for 



very small q, say q = A or 2 A. However, after carrying out the small-g limit for q m s (£), the 



discretization (|2lD of the integral in Eq. (|20D is harmless since the major contributions come 
from the intermediate wave vector domain k ~ 7. This reflects the fact that the sluggish 
dynamics of the tagged particle is ruled by the cage effect, i.e., by structure correlations on 
length scales of the interparticle distance. 

Second, there is a trivial relation between 5r 2 (t) and the velocity correlation function 
K s (t). Equations ( |I9]) and (|20| ) are equivalent to the MCT equation discussed earlier for 
K s (t) in the Laplace domain [[J . Following this route 6r 2 (t) has been evaluated previously 



for some representative packing fractions for the HSS ||26|| . We prefer to solve Eq. (|19D 
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directly, since thereby we need not worry about a careful handling of the strong small-s 
divergency exhibited by C [5r 2 (t)] (s). 



B. The diffusion localization transition 

In order to understand that the ideal liquid-glass transition implies a transition from 
particle localization to particle diffusion, let us remember that the Fourier back transform 
$ s (r, t) = (5 (r — [f(t) — t (0)])) of the tagged particle density correlator from the wave- 
vector to the displacement domain is the probability density for finding at time t the particle 
in a distance r from its starting position. In the glass this distribution levels off for long 
times at the normalized distribution f s (r). Hence the tagged particle is localized in the 
glass matrix. A characteristic localization length r s can be defined, for example, in terms of 
the small-g limit of the Lamb-Mossbauer factor: /* = 1 — (qr s ) 2 + O (g 4 ). This fixes the 
long-time asymptote of the mean-squared displacement for a > 

lim 5r 2 (t) = Qr 2 . (22) 

t— *oo 

Obviously, 5r 2 (t) = J § s (r } t)r 2 d 3 r } so that 6rf = J f s (r)r 2 d 3 r. Equation (fl9|) yields the 
general formula for r s in terms of sums over the glass form factors |§ 

r 2 s = I/Fmsd (/, / s ) • (23) 

In the liquid state the low-frequency behavior of the Laplace transform of all correlation 
functions and kernels is smooth ||33|| . Therefore, one derives from Eq. (|T9|) the well known 



formula for the linear divergency in time 

lim 5r 2 (t) jt = QD S . (24) 



Here D s is the tagged particle diffusivity. It is also called the long-time diffusivity in order to 
distinguish it from Dq, which determines the short-time asymptote lim t ^o 5r 2 (t) ft = QD^. 
From Eq. ( |I§D one readily derives for the ratio of D s and Dq 
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Note that this ratio is less than unity, which is an obvious manifestation of the cage effect. 

Figure ^| exhibits the evolution of the bifurcation dynamics as probed by Sr 2 (t) for the 
HSS. For very small t one observes short-time diffusion. With increasing density the short- 
time diffusion is suppressed due to the cage effect. The glass curves level off for long times 
at 6r 2 . With decreasing density the localization length increases up to some critical value 
r sc . The critical value r sc = 0.0746 fits nicely to the Lindemann melting criterion, as already 
noticed in Ref. 0. The liquid curves intersect the plateau r sc for times of order t a ; the 
time t a is indicated for the n = 9 curve by a dot in Fig. |fj. Then they leave the plateau 
according to von Schweidler's law in order to cross over to the diffusion limit 5r 2 (t) = 6D s t 
for times large compared to t' a . For the n = 9 curve, the latter asymptote is indicated 
by a dotted line in Fig. |6] and the time t' a is marked by a full square on the graph. The 
increase of Sr 2 (t) above 6r 2 is the a process of the mean-squared-displacement dynamics. 
It deals with the tagged particle's leaving of the cage. The initial part of this process, 
where 5r 2 (t) is close to r 2 c , is stretched over a large dynamical window if |e| is small. For 
times of order 20t' a , 5r 2 (t) has increased to about unity; then stochastic dynamics sets 
in and 5r 2 (t) follows the diffusion asymptote. Between the end of the regular short-time 
transient and the start of the a process, there is a mesoscopic window for another anomalous 
dynamics. It deals with the stretched approach of 5r 2 (t) towards the plateau 6r 2 . In this 
sense the bifurcation dynamics, i.e., the dynamics outside the transient, deals with a two- 
step relaxation process. Asymptotic expansions shall be used to describe the bifurcation 
scenario by analytical formulae, thereby providing an understanding of Fig. |6|. 



C. Asymptotic laws 
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1. P relaxation 

An equation connecting the Laplace transforms 5r 2 (s) and (s) of 5r 2 (t) and m^ ^ (t), 
respectively, within the structural relaxation regime t 3> to follows from Eq. flT3|): s<5r 2 (s) = 
6/ [sm' ' (s)j. By inserting the asymptotic expansions (|6]) and (|T4] ) into this equation one 
again obtains the general result for the /3 relaxation to next-to-leading order 

5r 2 (t) /6 = r 2 sc - h MSD G (t) - h MSD [h (t) + K MSD G (t) 2 + aK MSD ] ■ (26) 

Here r 2 c follows from Eq. ( p3|) with / and f s specialized to f c and f sc , respectively, and the 
other amplitudes read 

h MS D = r£ [FtisD (K n + Tl lSD (/ c , h')] , (27a) 



Kmsd = ~r~~ (/», /»') + ^msd / sc ) + 3~msd (A ^ s )] - (27b) 



Kmsd 



MSD l sc 

hMSD 



4 



"M5L> [ oas v / v / J r sc 

where we have exploited that the mode-coupling functional Tmsd (f c , f sc ) is linear in / and 
f s and introduced the short-hand notation like hK for (hK) = h p K p . 

In Fig. |7| the mean-squared displacement is compared to the asymptotic results in lead- 
ing and in next-to-leading order for a tagged particle in the HSS. The next -to-leading-order 
result explains the deviations from the leading-order one. The general trend of the devia- 
tions is again explained by the short-time corrections to the critical law and the long-time 
corrections to the von Schweidler law on the liquid side — or to the nonergodicity parameter 
6r 2 c on the glass side. The range of validity is smaller than for any of the q- vectors discussed 
in connection with Fig. [|, since the absolute value of Kmsd is larger than of all the Kf 
(i = 0, 1,2,3). The next-to-leading-order results are a considerable improvement over the 
leading ones: On the glass side, the second-order range of validity extends over almost three 
decades, while the first-order range — as defined at the end of Sect. [Ill C 1| — does not even 



exist. On the liquid side the next-to-leading order adds about a decade to the short-time 
side of the range of validity. It extends it beyond the boundaries of the figure on the long- 
time side. However this happens in part by accident because the exponent of the second 



term (t 2b , 2b = 1.16) in the von Schweidler series is close to the exact long-time exponent, 
given by the diffusion law (6D s t). 

2. a relaxation 

A description of the a process for the mean-squared-displacement dynamics can be 
derived from Eq. (|i~9l). Thereby one obtains the superposition principle 

Sr 2 (t) = 5f 2 (t/Q . (28) 

Here the master function 5r 2 is to be calculated from 

t 

J m^ c (t - t') <5f 2 (t') dt' = 6t, (29) 
o 

with the initial condition Sf 2 (t) = 6 [r 2 c + h^SD^] + O (t 2b )- The kernel is given by the 
mode-coupling functional at the critical point T^sd ano - ^ ^ ne density-fluctuation master 
functions: m^ c (t) = J= C MSD ($> (i) , <l s (t)V 

The various superposition principles imply coupling of the a-relaxation time scales or 
relaxation rates in the following sense 0]. Let us characterize the long-time decay of the 
variable A in the liquid by some time r^. This time diverges upon approaching the glass- 
transition: in the leading asymptotic limit for o — > 0— one finds ta = C^t'^. All times 
or rates are proportional to each other and follow a power law, specified by the exponent 
7: 1/t.a = |e| 7 , Eq. (0). The constants of proportionality Ca or Ta depend on the 
variable A and on the precise convention for the definition of ta- The scale coupling or a- 
scale universality is demonstrated in Fig. ^, where the diffusion coefficients of hard-sphere 
particles of diameters d s = 1.0 and d s = 0.6, and the a-scaling rate 1/t{* s of the density 
fluctuations for wave vector qi are compared to the asymptotic predictions. Although the 
asymptotic behavior is the same for all the quantities, D s for d s = 0.6 already deviates 
visibly from the asymptotic result for n = 4, while the other quantities start to deviate 
only for n = 3. This underlines again the nonuniversality of the deviations. Note that this 
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difference might not be due to asymptotic corrections in the sense discussed in this paper 
but could come from the mismatch of structural relaxation, governed by the MCT kernel 
m s q (t), and the transient dynamics, ruled by the times r q in Eq. (Q). 

In Fig. |5] the rescaled mean-squared displacement is compared to the a master function 
Sf 2 (t) for various packing fractions. The time scale t' a is taken from Eq. (^|). Note that this 
figure is a harder test of the asymptotics then just testing the time-density superposition 
principle, since it tests the scaling time t' a in addition to the shape of the master function. 



The leading correction to (|28|) is known to be of order |er| while the leading corrections 
to the factorization property are of relative order \/\o] ||18|| , which explains why a scaling 
works much better than (3 scaling. As we have already pointed out in Sect. |III C 2| the 



"quality" of a scaling is observable-independent — as can be seen by comparing Fig. || to 
Fig. |5] — even though first-order (3 scaling is much worse for 5r 2 (t) than for $^ (t). This 
is so because the dominant correction to the long-time (3 process is the corrections to the 
von Schweidler law (see Fig. |7|), which is absorbed into the short-time expansion of the a 
master function. The dashed and dotted curves in Fig. || show, how the deviations from the 
superposition principle for small rescaled times t are explained by Eq. ([18]). 

V. GAUSSIAN APPROXIMATION 

The small-g expansion of the density correlator which lead us to the equation 

of motion for the mean-squared displacement, can be viewed as the leading term of the 
Gaussian approximation [j36|-|38| 

(t) « e - q2Sr2(t)/6 . (30) 

The Gaussian approximation becomes exact for the short-time dynamics (either ballistic 
or — as in our case — diffusive motion) and for the long-time diffusion, and is known to work 



quite well for simple liquids such as Argon in its normal state [38]. So in the following we are 



primarily interested in the intermediate time window of the structural-relaxation regime. 
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In Fig. [TO some tagged-particle correlators as calculated for |e| = 10 7 / 3 are compared 



with the Gaussian approximation. The latter describes the behavior of the self-intermediate 
scattering function of the HSS reasonably for all considered wave numbers. Coming from 
small wave numbers, at which the Gaussian approximation is asymptotically exact, the 
deviations start to appear at the end of the /3-relaxation regime, i.e., in the von Schweidler 
region as exemplified for q = q . Going to still larger wave numbers like q2 and q^ the 
Gaussian approximation also deviates at the beginning of the /3-relaxation regime, i.e., in 
the critical region. 

The Gaussian approximation for the critical amplitudes can be found by substituting 
Eq. (|26|) into Eq. (|30|) and expanding the exponential 

/r^^-, (31a) 
h s w e~^V/>MSD, (31b) 
K\ « K MSD + l -q 2 h MSD , (31c) 



K\ « K MSD + - n - q 2 h MS D. (31d) 



1 

2 V / r 3 x 

Here K MSD has been introduced in analogy to Eq. (|TTD. These approximation results are 
compared to the various amplitudes in Fig. |TJ and 0. In both cases we find a good qualitative 
description for the amplitudes; the Gaussian approximation even describes the change of sign 
of K s and K*. The Gaussian results for f sc and h s are in better quantitative agreement 
with the exact results than the results for K* and K*, since the leading small-g corrections 
to ( ft lap and ( |31b| ) are of higher order than those for (|31c|) and ( 31d ). Obviously the results 



for d s = 1.0 are better than for the smaller particle with d s = 0.6. For the large particle 
they even give good quantitative descriptions for some of the amplitudes: practically no 
deviation can be found for the critical Lamb-Mossbauer factor f sc up to q = 15, and the 
critical amplitude h s is reproduced up to q = 5. 



These results also explain the deviations for the dynamics in Fig. [Lj]: Since the critical 
amplitude h s is overestimated by the Gaussian approximation, the dashed curves relax 
faster in the /3-relaxation region. On the glass side, a larger h s leads to a larger f s via 
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f s = f sc + h s \/ a I (1 — A). In the same way the corresponding figure for d s = 0.6 can be 
inferred from Fig 0. 



VI. NON GAUSSIAN PARAMETER 



A. The equation of motion 



A cumulant expansion can be used to systematically classify the deviations of the density 
correlator from its Gaussian approximation. The leading contribution to the deviations is 
proportional to the non-Gaussian parameter a 2 (t) . One gets 

^ (t) = e -^ 2 (')/6 j x + i aa {t) [q 2 5r 2 {t) /6 y + (g6) I (32) 

where a 2 = (3/5) [5r 4 (t) / 5r 2 (t) 2 ] - 1 with 5r 4 (t) = (|r (t) - f (0)| 4 ) HHH- Since the 
mean-quartic displacement 5r 4 (t) is proportional to the fourth Taylor coefficient in a small- 
q expansion of (t) = 1 — q 2 5r 2 (t) /3! + g 4 5r 4 (t) /5! + . . ., one can derive an equation of 
motion for a 2 from Eqs. (|9]) and (|19|) 

t t 
[1 + a 2 (t)] 5r 2 (tf + D S J m (0) (t - t') [1 + a 2 (t')} 5r 2 (t' f d£ = 6D S J [2 + m (2) (t - £')] 6r 2 (*') dt' . 



(33) 

Here a further mode-coupling functional for the kernel rrS 2 ^ (t) = J-'ngp (i) , $ s (£)) is 
introduced: 

^ (/• n = is? / nSM & + If ) <*■ < 34 > 



For a numerical calculation the integral (|34]) is expressed by a Riemann sum analogous to 
Eq. (|21|). The derivatives df^/dk and d 2 jf.jdk 2 are approximated through their numerical 
equivalents, e.g. df^/dk = (fk+A — fk-A) / (2A), except for g = 0.2 and g = 0.6, where 
they are calculated from a small-g extrapolation according to the Gaussian approximation. 
For the same reasons as discussed above in connection with Eqs. (0), (|20|) , and ( 21]) it is 
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preferable to solve the derived equation for a 2 (£) directly rather than to deduce the non- 
Gaussian parameter as the small-g limit from the numerical results for <M (t). We have 
nevertheless checked that the latter procedure can be followed provided one chooses a finer 
wave vector grid and a proper large-g regularization of the integral ([Top. 

B. The transition scenario 

To begin with, some comments concerning the general properties of a 2 (t) might be 
appropriate. Neither 5r 2 (t) nor 0:2 (t) are correlation functions and, therefore, the known 
general properties of positive definite functions like $ g (i) and (t) need not be valid 
for these quantities. With A(t) = \r(t) — r(0)| 2 denoting a positive observable, one gets 
Sr 2 (t) = (A (t)) and therefore Sr 2 (t) > 0. Since (t) = 1 - q 2 Sr 2 (t) /6 + O (g 4 ) one 
finds for completely monotone density correlators: — (—d/dt) 1 5r 2 (t) ^ for I = 1,2,... 
So within the structural relaxation window t > t , and for all t for our colloid model 
defined by Eqs. ([I]) and (|9|), the mean-squared displacement is an increasing function of 
time: (d/dt) 5r 2 (t) > 0. The inequality (A 2 ) > (A) 2 implies a 2 (t) > -2.5. Hence, the 
non-Gaussian parameter is bounded from below, but it can have either sign. At early 
times the density correlator exhibits Gaussian behavior, independent of whether Newtonian 
or Brownian dynamics is considered, and therefore a 2 (t — > 0) = 0. The long-time liquid 
correlators describe diffusion, which is a Gaussian process, and thus ct2 (t — > 00) = 0. From 
Eqs. (^) and ( |19D one finds for a < that ct2 (t — ► 00) = O (l/t). Hence, in the liquid state 
«2 (t) cannot be a monotonic function of time. Negative a 2 (t) means, that the probability 
for the particle to move very far is suppressed relative to the one expected for a random-walk 
process. Similarly, the probability for moving far is enhanced if a 2 (t) > 0. In the latter case 
the cage boundary is more fuzzy than in the former. 

Figure [11] exhibits the evolution of the non-Gaussian parameter 02 (t) upon crossing 
the bifurcation point for tagged spheres of diameters d s = 1.0 and d s = 0.6 in the HSS. 
For the glass state the curves exhibit arrest for late times a 2 (t — > 00) = /ngp = 1 + 



23 



Fngp (/, f s )- For the liquid state the a process manifests itself as a bump of the c^-versus- 
log 10 1 diagram starting with an increase above the plateau f^cpi reaching some maximum, 
and then decreasing to zero. The liquid curves for n > 10, i.e. for |e| < 10~ 3 , exhibit the 
superposition principle: a change of a causes a shift of the peak parallel to the logarithmic 
time axis without change of the shape. The first structural relaxation step deals with the 
approach towards f^ GP from below for t ^> to- The functions for |e| < 10~ 4 exhibit the 
two-step-relaxation scenario for t > 10 . 

For dense normal liquids a>2 is found to be positive in agreement with molecular dynamics 
results for liquid argon near its triple point. Our results for d s = 1.0 and n = 2 are of about 



the same magnitude as the value 0.13 found for this Lennard- Jones system fl38fl . Increasing 
the density towards the critical one implies an increase of a<i by about a factor three for both 
values of d s studied; but for the smaller spheres the non-Gaussian parameter is about three 
times larger than for the larger ones. One concludes that glassy dynamics in the liquid does 
not lead to dramatic changes of the magnitude of the non-Gaussian parameter, a finding 



also supported by Fig. ITU 



For long times, and for the smaller sphere for all times, a 2 is positive. This enhanced 
probability for the particle to move further is what one would expect as a result of the 
building of a backflow pattern in the liquid. However, for d s = 1.0 the negative plateau 
value f%GP implies, that for sufficiently large densities a dip to negative values for a 2 (t) 
appears. For d s = 0.6 the plateau is positive and the dip does not exist. The depth of the 
predicted dip, as opposed to the height of the maximum in the a regime, is not a structural 
relaxation phenomenon. Rather it is caused by the cross over from the short-time transient 
to the first structural relaxation step. The mode-coupling approximations are constructed 
to describe the long-time behavior due to the cage effect. The theory does not handle short- 
time collision effects correctly, and therefore the size of the dip predicted for d s = 1.0 might 
be an artifact due to the insufficiencies of the MCT in a regime, which it has not been made 
for. An analogous reservation applies, of course, to the quantitative details of the crossover 
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regimes displayed in Fig. |3| for (t) and in Fig. |6] for 5r 2 (t) . The MCT equations do not 
guarantee the validity of the inequality a 2 (t) > —2.5. Indeed for d s = 1.0 and if > <£> c + 0.01 
the inequality is violated and therefore the corresponding glass curves are not shown in Fig. 

The extension of the Gaussian approximation according to Eq. is illustrated by 



the dashed dotted lines in Fig. [10|. For the wave vectors q < q 2 , where the Gaussian 
approximation yields a good description of the correlators, the addition of the leading- 
cumulant correction proportional to a 2 (t) improves the fit seriously. However, for larger 
wave vectors the addition of the a^-term does not lead to improvements as is demonstrated 
for wave vector q 3 in Fig. |10|. As a 2 is considerably larger for d s = 0.6 than for d s = 1.0 the 
cumulant expansion (|32|) already breaks down for q near q . 

Figures ^|, [|, and |TT] have been shown in order to illustrate the theoretical essence of the 
MCT bifurcation. To avoid misleading conclusions from Fig. [TT] it might be adequate to 
remember the windows accessible by state-of-the-art experimental studies. In molecular- 



dynamics work 0,^] a variation of the diffusivity or other structural relaxation scales could 
be detected over a window of about four decades. This corresponds to the curves n < 7 in 
our figures. The dynamical window for structural relaxation explored by van Megen and 
Underwood H by photon-correlation spectroscopy is also about four decades wide. The 
size of the dynamical window accessible by the neutron-spin-echo instrument || is smaller 
than three decades. To focus on results for a 2 which might be relevant for the interpretation 
of experiments with the techniques available today, one should ignore the results for n > 8 



in Fig. |TT]. This restricted set of results does neither exhibit the superposition principle 
nor the two-step relaxation scenario. The leading-order asymptotic results for the MCT 



bifurcation dynamics do not describe the MCT results shown in Fig. [TT] for n < 8, i.e., for 
|e| > 0.001, not even qualitatively. But in the following it shall be shown that the results 
can be understood in terms of the next-to-leading asymptotic laws. 
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C. Asymptotic laws 



1. (3 relaxation 

In order to work out the dynamics in the /^-relaxation window we start from the Laplace 
transform of Eq. fl33|) 

C [(1 + a 2 (*)) 5r 2 (t) 2 } (s) = s5r 2 (sf {2 + sC [F NGP ($ (t) , <T (*))] (s)} , (35) 

where the Laplace transform of Eq. (|TJ5[) has been used. Inserting the asymptotic expansions 
(H), (0), and (EH) for $ g (i), $ s (t), and 5r 2 (t), respectively, one gets: 



« 2 (*) = Jngp + h NGP G (t) + h NGP \H (t) + K NGP G {tf + aiWj . (36) 

Here the amplitudes are given by 

h NGP = F C NGP (h, f sc ) + T C NGP (f c , h s ) , (37a) 
K NGP = [F C NGP (h, h s ) + T C NGP {hK, f sc ) + Fngp (f c , h s K s )] /h NGP 

+ hMSD_ ^ _ ^ ^2r 2 c h NGP - h M SD(l + /jvop)] /h NGP , (37b) 



sc 



K 



NCP 



+ T C NGP (hK, f sc ) + T C NGP (/ c , h s K> 



/h 



NCP 



df c NGP (f c , f s 
Cde 

+ hMSD_ [h MSD (i + f NGP ) _ 2r 2 c h NGP ] /h NGP . (37c) 

r sc 



The analytic results are compared with the numerical solutions for n = 9 in Fig. [12. 
Because of the large negative value of h NGP , the leading order qualitatively accounts for the 
steep rise of «2 in the /3-relaxation region. This is a major reason for the absence of a plateau 
even though the reduced packing fraction |e| = 0.001 is so small. The large size of Hngp 
also explains, why the long-time limit Jngp i n the glass is much more suppressed below the 
plateau value /nop, than shown in Fig. [7| for 5r 2 (t). There are two further peculiarities 
hidden in Eqs. ([H|), which render the discussion of Fig. [12] different from those considered so 
far for other examples. The a process for the density correlators and for the mean-squared 
displacement deals with monotonic functions. In those cases the leading corrections to the 
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von Schweidler law influence the details but not the general trend of the functions. The 
positive slope of (i) in the leading-order-/3-relaxation regime, however, has to change to 
a negative one in the late a-relaxation regime since the a.2 (t)-versus-t curve eventually has 
to decrease to zero. Therefore the «2 (t)-versus-log 10 1 diagram exhibits a bump. The next- 
to-leading asymptotic formula as opposed to the leading one, can reproduce at least 



qualitatively this new feature of the diagram as shown by the dot-dashed lines in Fig [12. 
The second difference concerns the remarkable parallel shift of the solution relative to the 
leading asymptote occurring for times near t~ where the liquid j3 correlator exhibits a zero, 
G {t~) — (t~ — 0.704t CT for the HSS). The correction is given by Iingp H (t~) + oKngp ■ 
For the previously discussed examples the two terms in the bracket partly cancelled each 
other while for a 2 (£) they add up. As a result this shift is about an order of magnitude 
larger than in the examples studied in Figs. |4] and [7| or in the examples considered in 



Ref. |18[ . This shift is responsible for the fact that for n = 9 the leading-order asymptotic 
result for ct2 (t) (d s = 1) (dashed lines in the lower panel of Fig. [12]) has a smaller range 
of validity than the corresponding result for the mean-square displacement (Fig. [1]) even 
though \K NGP \ < \K MSD \. Of course, this does not contradict the idea that the amplitude 
K asymptotically determines the range of validity of the leading-order result, because, 
asymptotically close to the critical point, the shift becomes irrelevant. Thus, the leading- 
order range of validity for a<i if) (d s = 1) will finally overtake that for Sr 2 (t) as a — > 0. 

One can substitute the expansions ( |26| ) and fl36|) into Eq. (^) in order to obtain the 
amplitudes in Eq. (|14]). The result extends the Gaussian approximation (|3TD for these 
quantities so that the first cumulant is taken care of. The result improves the Gaussian 
approximation for q < q$. However, for large wave vectors the extended approximation is 
worse than the Gaussian one, which is shown in Fig. |1] and |2| by the full lines. Therefore 
the described extension is not worthwhile. 
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2. a relaxation 



To explore the asymptotic behavior in the a-relaxation region, the non-Gaussian param- 
eter is plotted versus rescaled time t = t/t' a in Fig. [13]. Comparing the results for d s = 1.0 
with the corresponding ones for the mean-squared displacement in Fig. |^ one notices, that 
one has to choose a considerably smaller separation \a\ for a 2 it) than for 5r 2 (t) to find 
comparable agreement with the a-relaxation scaling law. For n = 5 and t — 1, 5r 2 (t) 
agrees with the superposition-principle asymptote within the accuracy of the drawing in 
Fig. |^, while the corresponding result for a 2 deviates seriously from the asymptote. The 
mean-squared displacement is dominated by the diffusion limit, i.e. Gaussian behavior for 



t > 20. Such contribution is absent in a 2 and therefore the n = 5 curve in Fig. [13| magnifies 
the small deviations from Gaussian dynamics, which for the n = 5 curve in Fig. ^ start to 
be visible only for t < 0.3. 

Figures [13] and || are seemingly in contradiction to the statement that a-scaling should 
work equally well for different quantities (see Sect. |III C 2|) . However, as mentioned above, 
the non-Gaussian parameter for the shown packing fractions is special in the sense, that 
the corrections for small t contain an important shift term, which stems from the term 



in brackets in the short-time expansion of the a-scaling correction (|T8p. In Fig. [13] this 
finding is illustrated for the curves with label n = 9. For even smaller \a\ this shift becomes 
irrelevant, thus resolving the paradox. 



VII. CONCLUSIONS 

In this paper the MCT results for the structural relaxation of the mean-squared displace- 
ment Sr 2 (t) and non-Gaussian parameter a 2 (£) were discussed for the hard-sphere system. 
More generally, we studied the generating function of these quantities, viz the incoherent 
intermediate scattering function or tagged-particle correlator (t). The work is motivated 
by the distinguished role of these quantities for the interpretation of spectroscopic data and 
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molecular-dynamics simulations of simple glass forming liquids. The essence of the MCT 
bifurcation scenario for the evolution of structural relaxation is contained in the leading 
asymptotic laws which deal with the two scaling laws reviewed in Sect. [H]. Most of the so 
far published tests of MCT focused on an assessment of these universal results. The outcome 
of this paper demonstrates in a drastic manner, that the range of validity of the mentioned 
asymptotic laws is not universal, rather it depends on the quantity considered. This range 
of validity can be determined by working out the laws for the leading corrections. Thereby 
formulae for a refined data analysis are obtained. There are quantities such as «2 {t), which 
are not described at all within the presently accessible dynamical window by the leading 
asymptotic result. In this case the laws including the leading-order corrections are necessary 
for a qualitative understanding. 

From Fig. ^ one infers that it is relatively easy to extract von Schweidler's law for 
the relaxation at the structure-factor peak position qi, while it is more difficult to identify 
this leading asymptotic law, Eq. (|4j), for the larger wave vector g 3 . The behavior of Q s q (t) 
for the intermediate wave-vector range qo-qs is quite similar to what was studied for the 



coherent intermediate scattering function in Ref. [IB]. The results, in particular those for the 
spectra, can therefore be inferred from that earlier work, using the amplitudes from Figs, 
[l] and ^. Inclusion of the calculated correction terms often extends the range of validity 
of the asymptotic formulae seriously. Therefore it is advisable to fit the initial part of the 
a process by the von Schweidler expansion including the P b term as was done recently by 



Sciortino et al. |25|] in their analysis of simulation data for supercooled water. These authors 



noticed in particular, that the leading correction term to Eq. ([|) was smallest near q\. This 
is in qualitative agreement with Fig. f|, even though there is no a priori reason to expect 
our hard-sphere-system results to be relevant for an explanation of the details of water. 

Fig. [7] implies that the 5r 2 (t)-versus-log 10 1 diagram exhibits only a small window for von 
Schweidler's law and that the window for the analytic description is extended significantly 
by inclusion of the t 2b correction term. It is also predicted that the critical decay law cannot 
be identified for reduced packing fractions |e| > 0.001. These findings are similar to what 
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was reported in Ref. |TT[ for the simulation results for a binary Lennard- Jones mixture. 



In agreement with photon-correlation-spectroscopy results for a hard-sphere colloid [22 



one infers from Fig. [H] and from the lower panel of Fig. [11], that the Gaussian approximation 
works very well for a hard-sphere system. This finding is not in conflict with the observation 
that 5r 2 (t) does not exhibit the critical decay law, while $> s (t) = exp [— q 2 5r 2 (t) /6] does 
show this leading asymptotic result, Eq. (§), for q ~ g 3 (see Fig. |). The strong t~ 2a 
corrections to Sr 2 (t) /6 ~ r 2 sc — Hmsd (to/t) a nearly cancel those coming from the expansion 
of exp [— q 2 Sr 2 (t) /6]; the relevant correction amplitude Kmsd for 5r 2 (t) is therefore larger 
than that for $3 (t) . 

The test of MCT against molecular-dynamics-simulation results by Kob and Andersen 



11,39,40|| gives strong support of the theory. This can be appreciated even more if one 



considers the comparisons of first-principle MCT calculations of the critical form factors, 



critical amplitudes, and exponent parameter with the data |41| . The authors also confirmed 
the power law singularity for the a-relaxation scales r oc \T — T c |~ 7 with the predicted 
connection between the exponent 7 and the von Schweidler exponent b, Eq. (|2|). However, 
they obtained the result for the diffusivity D oc \T — T c f with 7' < 7, in contradiction 
to the a-scale universality. This observation underlines in particular that the predicted a- 
scale coupling is not a triviality. Figure [8] demonstrates that the corrections to the leading 
asymptotic law for the diffusivity of a smaller sphere are larger than the ones for the a- 
relaxation scale for a representative intermediate wave vector q\. But our finding for the 
corrections to the hard-sphere-system asymptotics are too small to explain the reported 
results for the cited mixture. The MCT prediction for the universal a-scales for the HSS, 
T a oc [\ip — (p c \ /ip c ]~ 2 ' 6 oc 1/D, was confirmed by a recent analysis of the diffusivity D and of 



the tagged-particle-correlator a scale r a for a hard-sphere colloid by van Megen et al. |22 

The results for a 2 (t) in the upper panel of Fig. [11] are unusual if compared to the 
corresponding ones for representative density correlators $ q (t) or 3>g(t). For the curves 
with n < 8, dealing with the presently in experiments or simulations accessible parameters 
and windows, the two MCT scaling laws, Eqs. (j5|) and @, cannot be identified. Neither 
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is it possible to recognize the two time fractals, Eqs. (0) and (|J). These findings can 
be understood by analytical formulae only, if the next-to-leading asymptotic results are 



appreciated, as is shown in Fig. [L2[ However, these observations are not so surprising for 
the following reasons. The a process for the correlator $ 9 (t) deals with the monotonic decay 
from the plateau fi to zero. This phenomenon is described already within Maxwell's visco- 
elastic theory, dealing with exponential decay and Gaussian density fluctuations. Glassy 
dynamics deals with deviations from this phenomenological description. However, for the 
time range where Q s q (t) < fq C /2, which deals with a major part of the $^ (t)-versus-log 10 1 
graph, the phenomenological picture remains essentially valid. In this part the superposition 
principle holds already for n > 5, as can be inferred from the a-scaling analysis in Fig. |5| 
or from Ref . [18| . Also for «2 {t) the superposition principle works for the mentioned large 



times as shown in Fig. |13| But since there Gaussian dynamics is nearly valid, «2 (t) is 
small and the corresponding part of the a2-versus-log 10 1 diagram merely deals with a not 
so interesting feature of the figure. The diagram is dominated by the maximum and this is 
located near the end of the window of von Schweidler's law. Here the corrections are large 
compared to the leading terms, the maximum is only poorly approximated for n < 5. But 
Fig. |5] shows that for (t) the n = 5 results do not follow the a-scaling master curve for 
t/t' a < 0.2. This too can be inferred from the corresponding scaling plot shown as Fig. 17 in 
Ref. or from Fig. |9| for 5r 2 (t). In the latter case the deviations from scaling are judged 
relative to the underlying elementary background and Qr^ c , respectively, and therefore 
they do not appear as qualitative effects. The non-Gaussian parameter magnifies a small 
effect, and therefore its approximation by asymptotic formulae is judged differently than for 
the other functions. Indeed, for the HSS the predicted size a<i (t) < 0.3 is so small, that the 
effect disappears in the data noise of the measurement on colloids ||22|| . Let us emphasize that 



the value /ngp of the plateau was of importance for the preceding discussion. Therefore, 
it would be of great interest for an assessment of our results to measure /nop in the glass 
state. 



The upper panel of Fig. |TT] looks similar to what Kob and Andersen reported for their 
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results on a binary mixture |TT]. In particular we also find «2 to be larger for smaller 
particles; the maximum increases from 1.0 to 1.2 if d s is decreased from 0.6 to 0.5. The 
clustering of the a2~versus-log 10 1 graphs on the critical curve, labelled c in Fig. [U], and 
referred to as a scaling law in Ref. |fTI| , becomes better for d s = 0.5 than for the result shown 
for d s = 0.6. Let us emphasize, that the found results are not universal features of the MCT 
bifurcation dynamics. The apparent scaling gets disturbed if d s increases above 0.6 and it 



is predicted to be absent for d s = 1.0 as is shown in the lower panel of Fig. |TT 



It was shown by Kob et al. [fi2l, that the «2 peak is produced by clusters of particles 
which move faster than the ones in their neighborhood. Only as few as 5% of all particles are 
involved in the formation of these clusters. It would be very surprising if the MCT, which 
works with averaged quantities like $ g (t) and (£), could reproduce such subtlety of the 
microscopic dynamics. This holds even more so, if the speculations on a relation of the cited 



cluster dynamics with polymer dynamics jK| could be substantiated. On the other hand, the 



so far identified features of the fast clusters are not in an obvious contradiction to the idea 
that they are representative configurations building up backflow patterns. Moreover, the 
MCT equations for structural relaxation have originally been proposed and tested against 
experiments for a treatment of backflow phenomena Obviously it would be helpful 

to carry out a first-principle MCT calculation of «2 (£) for the cited binary mixture along 
the same lines as done in this paper for the hard-sphere system. Thereby one could clarify 
whether or not the qualitative agreement between the hard-sphere-system results for d s 



about 0.6 with the simulation results in Ref. [11| is just an accident. 
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FIGURES 

FIG. 1. Critical Lamb-Mossbauer factor /* c and amplitudes h s q , K^, K^, obtained from 
Eqs. (|l2|), (15), and (|l~7l), for a tagged hard-sphere particle of diameter d s = 1.0 immersed in 
a hard-sphere liquid. The solid lines are the Gaussian approximations, Eqs. ([H])- The arrows 
indicate the wave numbers go = 3.4, q\ = 7.0, q2 = 10.6, and q% = 17 A. The dotted line in the 
upper panel is one tenth of the structure factor at the critical point (compare Fig. 3 in Ref. ||i~8f). 
Here and in the following figures the diameter d of the host particles is chosen as the unit of length. 



FIG. 2. Same as Fig. |l] for a tagged particle of diameter d s = 0.6. 

FIG. 3. Self-intermediate scattering function 3>f(i), calculated from Eqs. @ and (|lT|). 
of the HSS {d s = 1.0) for wave number q\ = 7.0 at various reduced packing fractions 
e = (cp - tp c ) /ip c = ±10-™/ 3 , where n = 0, . . . , 14. The thick curve labelled by c, shows the 
dynamics at the critical packing fraction cp c = 0.515912.... The uppermost curve refers to the 
packing fraction ip = 0.6. The arrow marks the time to = 0.425. Here and in some of the following 
figures the full dot and square mark the times t a and t' a , respectively, for e = —0.001. On the glass 
side (e > 0) the curves for n = 13, 14 are omitted since they would be barely distinguishable from 
curve c. The n = 0, 1, 2 glass curves are omitted, since the Percus-Yevick approximation produces 
negative values for the pair distribution function for such high densities. Here and in the following 
figures the unit of time is chosen such that the short-time diffusion coefficient of the host particles 
reads D = 1/160. 
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FIG. 4. Upper panel: Rescaled self-intermediate scattering function 

(*) = {^q CO ~~ /g C ) /^o (solid lines) of the hard-sphere system (d s = 1.0) for four different 
wave numbers at reduced packing fraction e = (<p — <p c ) /<p c = ±10 _?1//3 with n = 9, showing col- 
lapse onto the (3 correlator G [dashed lines, Eq. (||)] in the /3-relaxation region. On the extreme 
left and right small numbers near each curve indicate the wave numbers qo, q%, q2 and q%, marked 
in Fig. |l|. Lower panel: Rescaled intermediate scattering function (ISF) and self-intermediate 
scattering function (SISF) for wave number q2 = 10.6 at n = 10, compared to the (3 correlator G. 
The glass curves (s > 0) are shifted downward by 1 to avoid overcrowding. The open symbols mark 



the boundaries of the /3-relaxation region as defined at the end of Sect. IIIC 1. Some symbols are 
missing because they are either outside the display range (the long-time part of the liquid curve 
for q\) or the range of validity is non-existent (glass curves for qo, q\ and q%). 

FIG. 5. a-scaling plot of the self-intermediate scattering function: $f (t), taken from Fig. 
||, for reduced packing fractions e = (ip — <p c ) /<p c = — 10~ n//3 , n = 1,3,5,7,9 versus rescaled time 
t = t/t' a . The thick solid line is the a master function ^ s q (Fj. The diamonds mark the early-time 



bound for the a-scaling regime as defined at the end of Sect. [II C 1. The dotted line indicates 



$g (t) — ah q Bii b , which includes the small-? correction according to Eq. (|l8|) . 

FIG. 6. Mean-square displacement 5r 2 (t), obtained from Eqs. (^) and (|2l|), for the 
hard-sphere system (d s = 1.0) at various packing fractions. The solid curves are labelled as 
in Fig. |!| The straight dotted line with unit slope indicates the long-time-diffusion asymptote 
(6D s t) for the n = 9 liquid curve. 
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FIG. 7. Mean-square-displacement curves 5r 2 (t) (full lines) taken for n = 9 from Fig. ^. 
The dashed and the dot-dashed lines are respectively the leading and the next-to-leading order 
of the asymptotic expansion ( |26| ) with r 2 c = 5.57 • 10 -3 , Yimsd = 0.0116, Kmsd = —1.23 and 
Kmsd = 3.33, calculated from Eqs. (|27|), The open symbols indicate the range of validity of 



the first (O) and second (O) order formulae as defined at the end of Sect. [II C 1. On the liquid 
side the second-order range of validity extends beyond the displayed range. On the glass side the 
first-order range of validity does not exist. The dotted straight line of slope unity for small times 
indicates the short-time diffusion asymptote (6-Dgi). 

FIG. 8. Self-diffusion coefficients D s for tagged particles of diameters d s = 1.0 and d s = 0.6 
and a-relaxation rate 1/t{* s versus reduced packing fraction e = (<p — ip c ) /(p c . The solid lines are 
the power-law asymptotes (r^ |e| 7 , 7 = 2.46). The a-scaling time t{* s is defined as the time when 
the tagged-particle density correlator <J?f (t), shown in Fig. ^ has relaxed to half of its critical 
plateau value ff c = 0.760. 

FIG. 9. a-scaling plot of the mean-square displacement: 5r 2 (t), taken from Fig. |(| for reduced 
packing fractions s = (ip — ip c ) /ip c = — 10~"/ 3 , n = 1,3, 5, 7, 9 versus rescaled time t = t/t' a . The 
thick solid line is the a master function 5r 2 (J) calculated from Eq. (|2^). The diamonds mark 



the early-time bound for the a-scaling regime as defined at the end of Sect. Ill C 1 . The dashed 
line shows the master function corrected according to the analogue of Eq. ( |i~8|) . Leaving aside the 
t-independent bracket term one gets the dotted curves. 

FIG. 10. Test of the Gaussian approximation ([}(]) (dashed lines) and of the cumulant expan- 
sion ( |32| ) (dot-dashed lines) for the self-intermediate scattering function of the HSS (solid lines) 
at packing fraction tp = (p c (l ± 10~ 7 / 3 ) for the wave numbers qo, q\, q2 and q%, introduced in Fig. 

a- 



FIG. 11. Non-Gaussian parameter ai (i), obtained from Eq. (|33|), for tagged particles of 
diameters d s = 1.0 and d s = 0.6. The labelling is done as in Fig. 0. 
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FIG. 12. /3-relaxation behavior of the non-Gaussian parameter Q2 (t) (full lines) taken 



for n = 9 from Fig. 11. The dashed and the dot-dashed lines are respectively the leading 
and the next-to-leading-order asymptotic expansion (|36|). The critical amplitudes (d s = 1.0: 
Ingp = -0.0445, h NGP = -1.12, K NGP = 0.543, K NGP = -2.98; d s = 0.6: f NGP = 0.891, 
hNGP = —1-37, Kn G p = 3.61, Kngp = —9.85) were calculated from Eqs. (|37[). The symbols 
indicate the range of validity of the leading (O) and the next-to-leading (Q) order as defined at 



the end of Sect. Ill C 1 — the first-order range of validity indicated in the lower panel is due to an 



"accidental" intersection of curves. 

FIG. 13. a-scaling plot of the non-Gaussian parameter: «2 (t), taken from Fig. O, for reduced 
packing fractions £ = (</? — Lp c ) /f c = ±10~ n//3 , n = 1,3,5,7,9 versus rescaled time t = t/t' a . The 
thick solid line is the a master function a.2 (t) ■ The dashed line shows the master function corrected 
according to the analogue of Eq. (|l~8|). Leaving aside the t-independent bracket term one gets the 
dotted curves. 
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